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Abstract 

As a function of the magnetic flux threading the object, the Little-Parks oscillation in the 
critical temperature of a large-radius, thin-walled superconducting ring or hollow cylinder has a 
period given by h/2e, due to the binding of electrons into Cooper pairs. On the other hand, 
the single-electron Aharonov-Bohm oscillation in the resistance or persistent current for a clean 
(i.e. ballistic) normal-state system having the same topological structure has a period given by 
h/e. A basic question is whether the Little-Parks oscillation changes its character, as the radius 
of the superconducting structure becomes smaller, and even comparable to the zero-temperature 
coherence length. We supplement a physical argument that the h/e oscillations should also be 
exhibited with a microscopic analysis of this regime, formulated in terms of the Gor'kov approach 
to BCS theory. We see that, as the radius of the ring is made smaller, an oscillation in the critical 
temperature of period h/e emerges, in addition to the usual Little-Parks /i/2e-period oscillation. 
We argue that in the clean limit there is a superconductor-normal transition at nonzero flux, as 
the ring radius becomes sufficiently small, and that the transition can be either continuous or 
discontinuous, depending on the radius and the external flux. In the dirty limit, we argue that the 
transition is rendered continuous, which results in continuous quantum phase transitions tuned by 
flux and radius. 
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I. INTRODUCTION 



The Little-Parks critical-temperature oscillations, with magnetic flux, of a large-radius 
thin- walled hollow cylindrical superconductor display a period h/2e l|, Q]. This oscillation 
period reflects the binding of electrons into Cooper pairs Qj. Recent experiments by Liu et 
al. 4j have probed the regime in which the diameter of the hollow cylinder is comparable to 
zero-temperature coherence of the superconductor. Their results confirmed the prediction 
by de Gennes j^J of the destruction of superconductivity in small rings for a certain regime 
of the external flux. Liu et al. raised the interesting and fundamental issues of what would 
happen if the circumference of the structure were to be smaller than the superconducting 
coherence length, as well as whether or not the Ginzburg-Landau approach would be valid 
in this regime. These issues make a microscopic treatment desirable. 

On the other hand, the single-electron Aharonov-Bohm oscillations in the resistance or 
persistent current in a clean metallic ring have period h/e [6|. This leads to a related 
fundamental issue: As the radius of the ring becomes smaller, how would the single-particle 
h/e period manifest itself in a Little-Parks type of experiment? Furthermore, does disorder 
affect the oscillation of critical temperature and the character of the transition between 
superconducting and normal states? And if so, how? These questions are not only of 
theoretical interest, but are also likely to be addressed experimentally, in view of recent 
progress in fabrication and experiments on small superconducting rings Q] • 

Recent work by Czajka et al. 9], involing the exact diagonalization of the Hubbard 
model for small numbers of sites and the numerical solution of the Bogoliubov-de Gennes 
(BdG) equation, shows that impurities can play an important role if they are located such 
that pinned density waves are in phase with one another: charge-density-wave (CDW) order 
would then be enhanced and superconducting order reduced. They also found that the mean- 
field results obtained via the BdG equations are consistent with the exact diagonalization 
results, even in the small systems they studied. Very recently, a numerical study by Loder 
et al. and analytical work by Juricic et al. 11 ] and by Barash [3| on clean <i-wave 
superconducting loops has shown /i/e-period oscillations in the supercurrent. 

Ginzburg-Landau (GL) theory is valid near the superconducting transition, as it is an 
expansion in powers of the superconducting order parameter. Therefore it may not be appli- 
cable at sufficiently low temperatures. Moreover, the GL approach cannot give a complete 
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account of multi-connected geometries of small size, as it is a description of the center-of- 
mass wavefunction of the Cooper pairs. For small rings, in addition to propagating around 
the circumference together, electrons in Cooper pairs can split apart and rejoin, and this 
process is not included in the GL description. In this Paper, we study the oscillations of the 
critical temperature of s-wave superconducting rings theoretically, via consideration of the 
microscopic BCS theory of superconductivity 13], analyzed using Gor'kov's approach fl^ . 
Our central purpose is to address the issues raised by Liu et al. 4j and mentioned above. 
Our focus is on the correction to the oscillations that is due to fmiteness of the radius of 
the superconducting structure. We consider both the clean and dirty regimes, and in the 
latter regime we shall ignore any tendency towards CDW ordering by assuming that the 
significant configurations of the impurities are sufficiently random that any pinned density 
waves are not in phase with one another. We shall see the emergence of an h/e oscillation 
in the critical temperature, as the radius becomes smaller, and we shall also see that the 
transition to the normal state can be either discontinuous or continuous in the clean limit 
but is always continuous in the dirty limit. 

The origin of the h/2e oscillations lies in the fact that a Cooper pair carries charge 2e, 
which when circumnavigating the ring acquires a phase e * 2e */ c?i , where $ is the magnetic flux, 
linking the ring. The period in flux is thus h/2e (with c conveniently set to unity [l5|). The 
emergence of h/e oscillations for small rings is due to the additional process in which electrons 
in a Cooper pair can, from time to time (so to speak), separate, propagate separately, and 
rejoin, the two trajectories having a nonzero winding number relative to one another. This 
process, which can only occur for ring sizes comparable to or smaller than the Cooper-pair 
size, induces an oscillation with period h/e. 

The organization of this Paper is as follows. In Sec. [Til we discuss the critical-temperature 
oscillations in the clean limit and in Sec. IHII we include the effect of disorder and discuss the 
dirty limit. We make some concluding remarks in Sec. IIVI In Appendix [A] we derive the 
effective one- dimensional Gor'kov equations by averaging over the cross-section (or thickness) 
of the ring. In Appendix [B] we provide supplementary details of the calculations that lead 
to the results in the main text. In Appendix [C] we provide two heuristic arguments for the 
emergence of the h/e period: (1) by examining Cooper's problem on a ring, and (2) by using 
an instanton approach. These two arguments lead to the physical picture described above, 
and complement the Gor'kov Green- function approach adopted in the main text. 
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FIG. 1: A ring with a magnetic flux threading through it. The radius of the ring is R and the 
thickness in the cross-section is d. 



II. CLEAN LIMIT 



A. Gor'kov equations 



We consider a ring (to be more precise, a torus) of radius R (and thickness d, namely, 
the diameter in the cross-section, which is smaller than both R and the Cooper-pair size £o) 
with a magnetic flux density B threading the ring parallel to the ring axis; see Fig. [TJ We 
can describe this field by the vector potential A(r) = (B/2)z x r. On the ring itself (the 
circumference of which is L = 2irR), the vector potential A is given by 9<&/L, where 9 is 
the azimuthal unit vector and $ is the total flux enclosed by the ring, i.e., <3> = J da ■ B(r). 

The normal and anomalous Green functions obey the Gor'kov equations [si, Q, 16] 



1 eA 

[ihw n - — (-iW - — ) 2 + fi}G{f, 7; uj n ) + A(f)F t (f, r*; u n ) = HS(r- r ), (la) 
1 e A 

[-ihw n - — (iW - — ) 2 + /i]Ft(r,r> n ) - A*(f)G(fy;u n ) = 0, (lb) 

where M is the electron mass, e is the electron charge, r and r 1 are three-dimensional 
coordinates in the ring, and the order parameter A is defined self-consistently via 

A ^) = ^E Ft ( r >~>«)' ( 2 ) 

UJn 

in which u n = 1-nTin + 1/2) are Matsubara frequencies, /3 = 1/ksT, T is the temperature, 
and V is the BCS pairing strength. 

We now invoke the narrowness of the ring to justify dropping all dependence on f and r* 
except that associated with the one- dimensional coordinates along the ring, x and x' . Owing 
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FIG. 2: Pairing of states [sj]. (a) Upper panel: the external flux is $ = 0; the pairing is between 
n \ + n>2 = 0. (b) Lower panel: $ = h/2e; the pairing is between n\ + n2 = —1. The pairing 
configuration changes from type (a) to type (b) at external flux = /i/4e; see Ref. 3]. 



to the physical periodicity of the ring, all functions of x and x' are periodic with period L 
(and the vector potential is a constant along the circumference of the ring). The Gor'kov 



equations then become one dimensional 



17|: 



+ %hu) n - (ihd x + ^) 2 + fj, G(x, x'\ uj ti ) + A(x)F j (x, x'; u n ) = hS(x - x'), (3a) 



itw n ~ ^(^^ _ + ^ F^(x,x';uj n ) - A*(x)G(x,x';uj n ) = 0. 



We can expand G and F in Fourier series, as follows: 
G(xi,X2',uj n ) = — g nitn2 (uj n 



(3b) 



x 1 ,x 2 ;uJn) 



ni,n 2 



'cjje i ^ nia:i+i ^' 12a;2 , 



(4a) 
(4b) 



where ri\ and n 2 are integers labeling single-particle states. Due to the translational (to be 



more precise, rotational) invariance of the system, we assume that G(x, x) has no depen- 
dence on x. This sets constraints on the nonzero Fourier components for G: U\ + n 2 = 0. 
Furthermore, we assume that A(x) = e l2nmx ^ L A , and hence F^(x,x) ~ e - l2nmx / L /\*. This 
sets constraints on the nonzero Fourier components for F*: U\ + n 2 = —m. The meaningof 
this is that the pairing occurs between the single-particle states ri\ and n 2 = —m — ri\ pj; 
see Fig. [21 

The Gor'kov equations can be expressed in terms of the Fourier components of G and F^ 
as follows: 

h 2 2i 
+ - 2MR2 {ni + m-4>) + /ij g ni+m - ni - m + A fl u _ ni _ m = H, (5a) 

fr 2 2l + 

~~ %fk ° n ~ OII/P2 ( Ul + 0) + ^ /ni,-m-m ~ A tfm+m -m-m = 0, (5b) 
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where = $/(— /ic/e) = $/(/ic/|e|). In the following we shall set ft = 1 , c = 1, and &b = 1, 
re sake of convenience. But we shall refer to the single-particle flux quantum ftc/e as 



for t 



ft/e 15]. These equations can be solved explicitly, yielding 



^ _ -iun ~ ^l{n 1 + (f)) 2 + fi ^ . 

[iuj n — VL(ni+m — (f)) 2 + //] [ — iu n — fl(rii + (f)) 2 + /x] + |A | 2 ' 

A* 

ft = f±5 f 6 b) 

Jni '" ni " m [iw n -fi(ni+m-0) a + ^][-ia; n -n(m+0) 2 + ^] +|Ao| 2 ' 1 ; 

where, for the sake of convenience, we have introduced Q = 1/2MR 2 (noting that we have 
set ft = 1 in h 2 /2MR 2 ). The self-consistency equation (j2J) then becomes 

VT 

We mainly discuss the case in which the chemical potential fi is kept fixed, e.g. by contact 
with a particle reservoir, or else we assume that the variation of /i with temperature and 
flux is sufficiently weak to be negligible near the superconducting transition. 



B. Critical temperature 

To solve for T c (<f>) we set A = in the self-consistency equation, thus obtaining 

VT 1 

L ^ ^ [iUn-VL^nx+m-cp) 2 + fi][-iu n -Vt(ni+<p) 2 + fi}' 
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It is important to note the underlying assumption that the transition from superconduct- 
ing to normal is associated with a vanishing order parameter, and hence is a continuous 
transition. The consistency of this assumption needs to be checked once we obtain the 
solution. We shall see, below, that for sufficiently small radii the assumption is not valid 
and the transition is actually associated with non-vanishing order parameter, and hence is 
discontinuous. 

Next, we make use of the root of the Poisson summation formula (i.e. the Dirac comb 
and its Fourier series), Y2 ni &z ~ n i) = ^2kez e l27Txh , to turn the summation over n\ in 
Eq. (jSJ) into one over the conjugate variable k: 



Instead of placing a cutoff on the energy, we follow Gor'kov [jjj] and place the cutoff (which, 
in practice, is the Debye frequency uo) on the Matsubara frequency. 

1. Large-radius limit 

For sufficiently large R we can ignore the correction terms associated with finite radius 
(i.e. k ^ 0), and thus we obtain an equation relating the critical temperature at nonzero 
flux to that at zero flux. In the limits that the Debye frequency is much smaller than the 
chemical potential (i.e. ujd/h *C 1), and that the chemical potential is much larger than the 
level spacing (i.e. \J 2MR 2 fi ^> 1), we obtain (see Appendix IB 21 for the derivation) 



where T c ° = T c (0) is the critical temperature at zero flux for the same radius, x m (<f)) = 
<j) — m/2 [with m being chosen to minimize (20 — m) 2 , i.e. the kinetic energy in the Ginzburg- 
Landau picture], and i/j(x) is the digamma function, 




(9) 




(10) 




(11) 



where 7 is the Euler constant 181 ] . 
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The zero-flux, zero-temperature order parameter Ao is related to the corresponding crit- 
ical temperature T c ° in the weak-coupling limit by Ao = rT c °, where r = n/e 1 ~ 1.76. We 
call the length £ = vf/kAq the Cooper-pair size (where vp is the Fermi velocity), so as 
to distinguish it from the zero-temperature Ginzburg-Landau coherence length £(0). There 
are two further relevant energy scales: the chemical potential /i, and the single-particle 
energy-level spacing Q. The ratio Qli/Aq can be estimated as 

Sf-(^)(^)/(S , -@)'- 

which is set by on the ratio of the ring radius to the Cooper-pair size. This motivates us to 
define a measure p of the ratio of the radius to the Cooper-pair size, via 

7T \ 2 _ Q jl 



(13) 



>y a 2 , 

Then, defining t = T c /T c °, we can re-write the equation for the critical temperature as 

lni ^(0_ Re ,Q_-^Y (14) 



An equation of this form was studied by Sarma [19j and by Maki and Tsuneto [20|, both 
in the context of the effect of a magnetic exchange field on superconductivity. The cor- 
respondence is that the role of the exchange field normalized to the zero-temperature gap 
(i.e. liH/hAq) is, in the present setting, played by the combination of the normalized inverse 
ring radius and the flux (i.e. x m (4>)/p). In Refs. Q,Q it was found that for large enough 
exchange field (here, small enough ring radius) Eq. (fl4"l) has multiple solutions and that, 
moreover, the correct interpretation is that the transition between the normal and super- 
conducting states becomes discontinuous (beyond certain value of the exchange field). By 
contrast, for small enough exchange field (here, lar ge e nough ring radius) the transition is 
continuous. By borrowing the results of Refs. 19|, |20(], as summarized in Fig. [3j we have 
that the threshold at which the transition changes character between continuous and dis- 
continuous occurs at \x m (4>)\/p ~ 0.6/V (i.e. point B in Fig. [3]). Furthermore, one has that 
for extremely small rings, such that \x m ((j))\/p > \/yj2it (i.e. point D), the system never 
becomes superconducting 2l|. Of course, these estimates would need to be modified if the 
finiteness of the radius were to be taken into account, but we expect that the qualitative 
separation into discontinuous- and continuous-transition regimes would still hold. 
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FIG. 3: Phase diagram: reduced critical temperature t(4>,R) = T c ((f), R)/T c (0, R) vs. flux and 
radius v = 7r\(f> — m/2\(£o/R) in the clean limit and in the limit that the finiteness corrections are 



19. 



20J|, the 



ignored (i.e. bulk limit). In the case of the exchange-field effect, discussed in Refs. 
horizontal axis becomes v = fi^H/Ao, where /zg is Bohr magneton and H is the exchange field. 
The curve ABC is the solution to Eq. f)14[) : ABD is the curve representing the true equilibrium 
transition line. From A to B the transition is continuous, whereas from B to D the transition 



is discontinuous, a construction first made in Refs. 



19, 



20 ] in the context of the exchange-field 



effect. The curve BE represents the metastability limit for "superheating", whereas the curve BC 



represents the metastability limit for "supercooling" 



For large-radius rings (i.e. p ^> 1), we can expand the second digamma function in 
Eq. (I14p to second order and the logarithm to first order, thus obtaining 

_ (1 _ i)K ^Q)(^) 2 ,_, 41 (^) 2 . (15) 

Hence, we see that the fractional reduction in the critical temperature due to the flux is 
given by 

(1 _ () „ SM (^) 2 „ Ml f^)' = Ml (£)' (2, - m f. ,16, 

The integer m must to be chosen such that (20 — m) 2 is minimum, so as to obtain the most 
stable solution; hence, we recover the standard Little-Parks oscillation result, for which the 
period is h/2e. Equation (flEI) can be re-expressed as 



t AR . , 

min \2<p — m\. (17) 



8.41 r^ rnez 

The l.h.s. can be much smaller than unity if R is much smaller than £o> but the value of the 
r.h.s. can range from to 1/2, depending on the value of external flux. This means that 

9 
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FIG. 4: Critical temperature (normalized to its zero-flux value) t vs. flux eft = $/(/i/e) for R = 5£o 
(upper) and R = 1.5 £o (lower) in the clean limit. The plots are made with cos (27m) = 1 and with 
only the k = 1 term in Eq. (jl9[) retained. 



there can exist a range of fluxes for which no solution of t exists. This reflects the fact that 
superconductivity is destroyed over certain ranges of flux 5|. To make connection with the 
result by de Gennes {5], let us multiply Eq. (fT7|) by 2ty and take the cosine of both sides. 
We recover the de Gennes result for the transition temperature (for the case in which the 



length of the side arm in Ref. 



5H is set to zero 
R 



cos 2tx 



cos ( 2ix 



(18) 



^(t) J V h/2e, 

where ~ 0.74r^o/\/l — i is the temperature-dependent Ginzburg-Landau coherence 
length in the clean limit. We note that, however, the Ginzburg-Landau approach is, strictly 
speaking, valid only near t ~ 1. Furthermore, the existence of discontinuous transitions is 
beyond the reach of the Ginzburg-Landau approach. 
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FIG. 5: Reduced critical temperature t{4>) = T C (0)/T C (O) (normalized to its zero-flux value) vs. flux 
4> = 3>/(h/e) for R = 1.25 with cos (2/c7rn) = 1 in the clean limit. The higher T c branches (blue) 
are calculated by retaining only the first two terms in Eq. (|19p . whereas the lower T c branches 
(red) are calculated by retaining as many as ten such terms. The plot shows multiple solutions to 
Eq. (|19p . and signifies the emergence of discontinuous transitions. The upper and lower branches 
of the solutions are expected to merge at certain values of </>; the appearance of a gap between 
them is an artifact of our considering only a finite number of values of <f>. 



2. 



Finite-radius correction 



What is the correction to T c that arises from the finiteness of the radius? By taking into 
account this correction we arrive at the following equation obeyed by the critical temperature 
(see Appendix IB 31 for the derivation): 



lnt(0) 



T_ zx m (0)r 

1 2 2pt((j)) 



-Re 



- 4 {cos(27rA;rT) f c (0, p, 1) 

k=l 

" fc "e-^/ c (x m (0),p,t(0)) cos 2^0] } , (19) 



where we recall that t(<ft) = T C (0)/T C (O), T is the BCS constant r = ir/e 1 ~ 1.76, we have 



defined n = y/2MR 2 fi, which is related to the number o 
is defined via the hypergeometric function 2-^1(0, b; c; z) 18l |: 

, / ,js j\ „ /I ix m ((f))T 3 ix m ((j))T _i E k £ t 

fc{x m {(P),P,t) = 2^1 " 



Cooper pairs, and the function f c 



'(l- ' Xm(0)r ) (20) 



2pt 2 2pt )' v pt 

In Figs. H] and [5] we show the flux dependence of the critical temperature for two particular 
ring radii. For the larger radius case, the value of t at $ = h/2e is essentially the same as 
that at $ = h/e and 0; for the smaller radius case and for cos(27rn) = 1 (i.e. all the pair 
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states are occupied at and below Fermi level), the amplitude at $ = h/2e is reduced, relative 
to that at $ = h/e and 0. Thus, as the radius becomes small we clearly see the emergence 
of the single-particle flux quantum period h/e. Also worth noticing is the occurrence of a 
second solution (with lower value) for the critical temperature at sufficiently small radii, 
shown in Fig. Compared to the higher T c solutions, for which our approximation of the 
T c equation (Tl9|) . by keeping only the first two correction terms yields rather precise values, 
the evaluation of T c at these lower values of temperatures requires more terms (e.g. four or 
more) to be included in order to maintain the same accuracy. 

For the case of large rings (i.e. p ^> 1), to determine the leading corrections it is adequate 
to retain only the k = 1 term and to set the correction ~ ix m (4>)T / pt to zero, when compared 
to values of order unity, in the second f c function. Thus, we arrive at the formula 

lnt(0)« -^^(h-i^f) -4cos(27rrl)e- 2 ^/ r /(0,p,l)[l - cos2tt0]. (21) 

The change in the reduced critical temperature t((f) is then approximately given by 

1 -t{(j>) « 8.41 (Jpj (20-m) 2 + 4cos(27rrl)tanh" 1 (e"^)(l-cos27T0), (22) 

where we have used the fact that 2-Fi(l/2, 1, 3/2, y 2 ) = tanh _1 (y)/y. Again, the integer m is 
to be chosen to minimize (20 — m) 2 . We see that, in addition to the parabolic dependence on 
4>, there is a sinusoidal correction of period h/e. This is the emergence of the single-particle 
flux dependence. We also note that this correction is not universal, in that it depends 
sensitively on the value of /i [and, moreover, the form of cos (27m) results from the simple 
quadratic single-particle spectrum]. It can happen that the correction due to the finiteness 
of the radius actually increases the critical temperature, i.e. when cos (27m) < 0. 



As we have argued using the results of Sarma |19|] and Maki and Tsuneto [20], the 
occurrence of multiple solutions in Eq. (Ti"4"l) for T c (0), for certain ranges of \x m (<j)) \/p, leads to 
a change from a continuous to a discontinuous superconducting-to-normal phase transition. 
Even with the corrections to T c due to the finite-radius effect, as the radius of the ring 
decreases (to a value comparable to the coherence length), we observe that Eq. (TIT)]) still 
possesses multiple solutions for T c (0) for certain ranges of \x m (<f))\/p ; see Fig. [5j This 
implies that somewhere in these ranges (of \x m (<j))\/ p) there exists a change from continuous 
(at larger radius) to discontinuous (at smaller radius) superconducting-to-normal transition. 
This is shown schematically in Fig. [61 If the radius is sufficiently large, the transition 

12 



t 



qT.s 



s 




\ 



\ 



J.2 



i 10 

-0.5-0.25 
t 



0.25 0.5 




i i \ I I 
-0.5-0.25 



0.25 0.5 



\ 




FIG. 6: Schematic depiction of the reduced critical temperature t(4>) = T C (0)/T C (O) (normalized 
to its zero-flux value) vs. flux <f> = <3?/(/i/e) for small radius rings in the clean limit. As illustrated 
in Fig. [5l for sufficiently small radii, there are multiple solutions for T c , as exemplified there 
by the higher (blue) branches and the lower (red) branches. Near (f> = and ±1/2, the upper 
branches are the equilibrium phase boundary and the transition to normal state is continuous. 
Away from these regions, the globally stable equilibrium phases must be sought by free-energy 
consideration, and the correponding phase boundaries are indicated schematically by the solid 
lines and represent discontinuous transitions. Upper panel: For all flux values, there exists a 
superconducting state. Lower panel: For smaller radius, it can happen that there are flux values 
for which no superconducting state exists. 

is always continuous. If the radius is sufficiently small, the the curve representing the 
discontinuous transition in one "dome" can intersect with that of the nearby dome (above the 
"void" region where no solution for T c of Eq. (fT9j) exists). If this void region is large, the curve 
of the discontinuous transition can go to t = at certain value of <fi without intersecting that 
from the nearby dome. However, calculations of T c (0) alone cannot determine the precise 
location of the change from continuous to discontinuous; considerations of free energies are 
necessary to settle this issue. 
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FIG. 7: Critical temperature (normalized to its zero-flux value) for the disordered regime vs. flux 
cf) = 3>/(/i/e) for R = 0.2£o and the mean-free path l e = 0.2£o (solid blue line) and for the same 
radius but l e <C R (dashed red line). The plots are made with cos (27r/cn) = 1 and (— l) k for the 
upper and lower panels, respectively. Moreover, only the k = 0, 1 and 2 terms in Eq. f|33[) have 
been retained, as the convergence is rather good. 

III. DISORDERED REGIME 

How does disorder, in the form of potential scattering from fixed impurities, alter the 
physical picture that we have obtained so far? Does it affect the critical temperature? Does 
it change the order of the superconducting-to-normal transition? We now address these 
issues, assuming that the configuration of the impurities is sufficiently random that the 
tendency towards CDW formation is not enhanced. 



A. Expansion of anomalous Green function; disorder average 

For the most part, we shall consider impurities that produce scalar potential scattering; 
mutatis mutandis, the effects of exchange and spin-orbit scattering can be straightforwardly 
included. To begin with, we include potential scattering via the V(x) term in the Gor'kov 
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equations |14J, which now read (here we restore the constants fi and c) 



ihLO n - —(ihd x + —) -V(x) + /j, G(x,x';u n ) + A(x)F 



M(x-x')(23a) 



2M 

i^n~ ^{iftd x - —) 2 -V(x) + [i F\x,x';u n ) - A*(x)G(x,x';uj n ) = 0, (23b) 
where ^(x) represents the potential from static impurities, i.e., 



V{x) = ^u(x- x a ), 



(24) 



and x a indicates the spatial location of a th impurity. 

Following Gor'kov 141 ]. we introduce a Green function G°(x, x';u n ) that satisfies 

- ^^(^^ + ~~ + G Q {x,x'\ w„) = M(x - a;'). 
We can then express exactly in terms of A* and G° as 



F\x u x 2 ]w n ) 



h 



dx G°(x, x\\ u>n) A*(x) G(x, x 2 ; —uj u )- 



(25) 



(26) 



Because we are only concerned with solving for T c , near the transition, it is sufficient to keep 
G to zeroth order in A and replace G in Eq. fT2"6"j) by G°; thus we have 



F j (xi,x 2 ; uj n ) « ^ y dxG°(x, xi,u n ) A*(x) G°(x,x 2 ; -u n ). 



(27) 



We now consider the self-consistency equation ([2]), and average over the quenched disorder 
associated with the locations of the impurities: 

where 7TT indicates disorder averaging. We note that, as explained by Gor'kov jjjj], the 
Green function G° oscillates on a much smaller length scale than A* and, hence, the disorder 
average of A* can be factorized. For convenience, we use G to denote the disorder average 
of G°, i.e., G(x, r;u n ) = G°(x,r;uJ n ), the translational invariance of which is restored, viz., 



G(x,r;u} T 



t—r-L(x-r) 



(29) 



To calculate the disorder average of the product of two Green functions, we introduce 

n 

the kernel K [14J, defined via 



G%x,r;u n )G%x,r';-u n ) = ^ ^ K t 



. Zirn-\ , \ . Z7T710 , i\ 



(30) 
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If we retain only the ladder diagrams (i.e. ignoring the crossed diagrams 22f] ) , we arrive at 
the result 



K LJ {n u n 2 ) = G w (ni) G-^i^) [l + n imp ^ \u(q)\ 2 K w (m - q,n 2 + q)], (31) 

i 

where we have simplified the notation by dropping the subscript nonw and moving uj from 
an argument to a subscript, and n imp is the impurity concentration. Assuming that the 
order parameter retains the form A*(r) = Aq e ~ l2nmr / L ^ the self-consistency equation can be 
reduced to 

1 = 777 5 n 1 +n 2 ,mK UJ (n 1 ,n2). (32) 

ui ni,ri2 



B. Critical temperature 



We leave the detailed calculation of the kernel to Appendix IB 41 and simply quote here the 
resulting equation (IB48I) for the critical temperature in the disordered regime, i.e. Tq T c <C 1: 



hit 



j + 2_^^ e cos{2irkn) x 



k=l 



2Tvk R 

r 



^008(2^0)^(^(0), R,t{cf>)) -e —^f d (0,R,l) 



(33) 



where l e = vpTQ is the elastic mean free path, To is the elastic scattering time, we recall that 
T = 7r/e 7 , and 



f d (xM),R,t(<l>))= 2 F 1 



1,- 



r So 



■(l + HE^L). (34) 



ti? 2 ' '2 tR 2 

We remark that the argument in the second digamma function in Eq. (133]) is real, in contrast 
with the clean case, for which it is complex; see Eq. (TlOl) . A consequence of this is that there 
is no longer a doublet of solutions for T c . Moreover, the resulting single solution is consistent 
with the assumption that the order parameter becomes vanishingly small as the temperature 
approaches its critical value from the superconducting side. Therefore, the transition to the 
superconducting state is continuous. 

To explore the consequences of Eq. (1331) . we first examine the large-radius (i.e. bulk) limit, 
and show that we recover the de Gennes results for the case of rings. We then proceed to 
compare how the finiteness of the radius affects the oscillations of T c (0). 



1(3 



0.5 1 1.5 2 



FIG. 8: Phase boundaries between normal state (to the left of the curves) and the superconducting 
state in the temperature (t((p,R) = T c ((p, R)/T c {(p = 0,R), normalized to its zero-flux value) 
vs. radius (R, in unit of £o) plane. Upper panel: (f) = &/{h/e) = 0,0.1,0.2,0.25 (from the top 
down). The solid blue boundaries take into account of the finite-radius corrections, whereas the 
dashed red boundaries are solutions to the bulk equation, Eq. (i35j) . Lower panel: 4> = &/(h/ e ) = 
0.5, 0.40, 0.30, 0.25 (from the top down). In all plots, we have chosen l e = 0.2£o and cos (27r/cn) = 1, 
and have kept terms up to k = 2 in Eq. (|33p . Note the significant deviations associated with the 
finite-radius corrections near <p = 0.5. 

1. Large-radius limit 

Ignoring correction due to the finiteness of the radius, we obtain the following equation 
for T c of the bulk superconductor, i.e., 



as we should. This is in a gree ment with the results of de Gennes' j^J for rings and those 



obtained by Lopatin et al. 24| for hollow cylinders. To see that the de Gennes results are 
recovered, we note that the critical value of Yl P ^x^(d)) / , beyond which no superconduct- 



ing solution for any T > exists, is given by T /An 
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251 ] . For a fixed radius, this defines the 



critical flux $ c , which defines the boundary between normal and superconducting states. 
We then determine that the critical flux <p c = <& c /(h/e) satisfies 

— ^^=- = 27r min 1 20 c — m | , (36) 

which gives the critical flux for a given radius, or vice versa, via 



cos 



R ^ oos(^). (37) 



When the flux dependence of the critical temperature is weak, as in the large-radius limit, 
t = 1 — T c /T® is close to unity and Tl^QX^cj)) / R 2 <C 1, so one can expand the logarithm in 
(1 — t) and the second digamma function around 1/2 in Eq. fl35l) to obtain 

(38) 

As we did for Eq. (fIBl) in the clean limit, we can express this equation as 

cos ( 2vr-^- ] » cos f 2tt— — J , (39) 



where £(£) ~ 0.84rv£oII/\/l — t is the temperature-dependent Ginzburg-Landau coherence 
length in the dirty limit. Thus, we recover the de Gennes results but via a microscopic 
calculation. We note that, strictly speaking, the Ginzburg-Landau approach is only valid 
near t ~ 1 (see also discussion at the end of Sec. Ill B II) . 



2. Finite-radius regime 

Strictly speaking, the de Gennes results (13"9"j) are only valid in the large-radius limit, as is 
clearly seen from our microscopic derivation. In the regime in which the radius of the ring 
is comparable to the zero-temperature coherence length or the Cooper-pair size, we should 
take into account the effect of finite radius and use Eq. (|33|) instead of Eq. (135|) . In addition 
to the lengthscale defined by the Cooper-pair size £o, the behavior of T c also depends on 
another length scale, viz., the mean-free path l e . The correction terms (i.e. the terms in the 
summation) in Eq. (I3"3"|) decay exponentially with R/l e , so that the series converges rather 
rapidly, even for R ~ l e or slightly smaller. In Fig. [7] we contrast the predictions of T c in the 
large-radius limit, Eq. (135|) . to those that include the finiteness corrections, Eq. (1331) . for a 
small radius (R ~ 0.2£ , h)- Figure [7] also shows that the T c (0) oscillation has a component 
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of period h/e, in addition to the usual Little-Parks component of period h/2e. In Fig. [8] we 
plot the critical temperature vs. radius for various flux values. We see that the deviation 
from the bulk result is mroe significant near = 1/2, whereas the bulk result shows almost 
no deviation in the range <p = — 0.25. 

The deviation of the critical temperature of small-radius rings from that of large-radius 
rings is, however, not universal, as it depends on both the magnitude and sign of cos(27m) = 
cos(2tt ^2MR 2 fi). We have seen that when the sign of cos (2n a/2 MR 2 ji) is positive, the 
corrections from the finiteness of the radius can cause a reduction in the critical temperature 
near the flux value $ = h/2e, compared to the bulk case. This corresponds to the case in 
which, at zero flux, all the electrons are paired up. However, when the sign of the cosine is 
negative, the critical temperature at $ = h/2e can actually exceed its value at zero flux; this 
corresponds to the case in which, at $ = h/2e, all electrons are paired up, thus the system 
is more stable there than at $ = 0. If the cosine term should happen to vanish accidentally, 
the oscillation at period h/e would disappear. 

When the radius is not small the change in T c (<j)), relative to T c (0), is small, and we 
therefore need only retain terms in Eq. ( l33i) to order k = 1 and also may replace t on the 
r.h.s. by 1, as we did to obtain Eq. fl38|) . Thus we obtain 



Hence, an oscillation of period h/e is clearly seen to emerge, and whether the critical tem- 
perature reduces or increases (e.g. at $ = h/2e), compared to the bulk Little-Parks value, 
is seen to depend on the sign of the cosine term. For the case where all electrons are paired 
up at zero flux, the critical temperature is lower at $ = h/2e than at $ = 0. 

IV. CONCLUDING REMARKS 

We have considered the oscillations in the critical temperature of a superconducting ring 
of finite radius in the presence of a threading magnetic flux. We have found that, as the 
radius of the ring is (parametrically) reduced, an oscillation in the critical temperature of 
period of h/e emerges, in addition to the usual Little-Parks dependence (the period of which 
is h/2e). Our results provide corrections, due to the finiteness of the ring radius, to the results 
that de Gennes obtained for a flux-threaded ring [sj]. We have argued that in the clean limit 



(1-f) 



7T 



n e g 

2 R 2 
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there is a superconductor-normal transition, as the ring radius becomes sufficiently small at 
nonzero flux, and that the transition can be either continuous or discontinuous, depending 
on the radius and/or flux. In the disordered regime, we have argued that the transition is 
rendered continuous, which results in a quantum critical point tuned by flux and radius. 

One may wonder how the system behaves as it goes from clean to dirty limit. At which 
point does the existence of multiple solutions in the critical temperature disappear? By 
analyzing Eq. (IB39I) . we obtain that, ignoring the finiteness corrections, double solutions 
disappear when the disorder is such that l e /£o ^ 1-73. 

One question we should also address is the thickness d of the ring cross-section. It causes 

an orbital pair-breaking effect. For the purpose of estimating this we can use the result from 

a wire with same thickness in a field perpendicular to the wire axis (for the calculation, 
fl 

see, e.g. Ref. |2|). The fractional decrease of critical temperature, when it is small, can be 
estimated to be (not including the oscillation by flux) 

which causes a quadratic decrease in the critical temperature on top of the oscillations 
discussed in the present paper. If we take the same values shown in Fig. i.e., R ~ l e ~ 
0.2£o, for d ~ 0.2R the decrease is about 3.6% at <fi = 1/2 and about 14.4% at <p = 1. 
The largest decrease in the critical temperature shown in Fig. [7] is about 20%, and the h/e 
component is not swamped by the pair-breaking effect and can still be observed. 

Although we have obtained the emergence of the h/e oscillation and its amplitude by 
the microscopic calculations, the physics behind these effects can be illustrated via heuristic 
arguments associated with the corresponding Cooper problem on a ring, together with a 
path-integral based instanton tunneling approach. We briefly discuss these points in Ap- 
pendix o 



The non-universal factor cos(27rn) = cos {2-k^/2MR 2 ^ determines whether the finiteness 
of the ring radius leads to a decrease or an increase in the critical temperature near flux 
$ = h/2e. Although this factor is model dependent (i.e. dependent on the form of the 
single-particle spectrum), for the quadratic spectrum we consider here, \/2MR 2 fi roughly 
counts the number of electrons divided by four, as there are two spin species and positive 
and negative (angular) momenta. If we restrict ourselves to the case in which all electrons 
are paired then n = yj2MR?ii ~ (N pa , ir — l)/2, because for /i — two electrons can 
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still occupy the n = state. When all the levels at \i and below are filled at zero flux, 
cos (27r v /2M J R 2 /i) « 1. When there is one pair fewer or more (or equivalently, when all 
pairs are occupied at $ = h/2e), cos (2vr ^2MlPu) m -1. This seems to result in an 

261 ] but from the number of Cooper pairs. 



even-odd effect, not from the number of electrons 



However, whether this even-odd effect holds in general requires further investigation. 



Notes added in proofs. 1. We would like to point out that the issue o 
supercurrents for s-wave rings was also studied by Zhu and Wang in Ref. 



flux- dependent 
. In addition to 



Refs. 



10 



in Ref. [29 1 



11 



121 ]. there is a more recent work on the same issue in d-wave superconductors 



2. Throughout our paper, we have essentially assumed that the switching of pairing 
configuration (see Fig. 2) occurs at flux values being an odd integer multiple of h/4e near 
superconducting-normal transition temperatures. This is consistent with our numerical re- 
sults that by allowing the switching to be varied, the largest possible critical temperature 
is obtained when the above assumption is obeyed. In a very recent paper by Vakaryuk in 
Ref. 30J], the author concludes that the switching can be flux dependent at zero tempera- 
ture. If the two results are to be consistent, we are led to the conclusion that the switching 
is temperature dependent. This requires further investigation. 
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APPENDIX A: AVERAGED GOR'KOV EQUATIONS OVER RADIAL DIREC- 
TION 

Consider a ring on a plane with inner and outer radii a and A, respectively. Assume 
that the thickness (A — a) is much smaller than the mean radius (A + a)/2 and the zero- 
temperature coherence length £(0). We reduce the two dimensional Gor'kov equations into 
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effective one-dimensional equations by defining the averaged Green function, 

= ^^J A dp P J dp'p'G(p,e,p',9';u), (Al) 

and similarly for F\ Under this averaging, the two-dimensional delta function becomes 
one-dimensional, 

pA pA 

/ dpp dp'p'5(r-f') = 5(9-9'). (A2) 

J a J a 



A 2 - a? 



In the Laplacian operator, there is a term --^-p-g-G(p, p'). When performing the above 
average, we get for this term 

f A f A 1 d d , 2 f A , d p= a 

/ dpp dp'p'——p—G(p,9,p',0';uj) = — ^ dp'p'p—G , (A3) 

J a J a P dp dp A 2 - a 2 J a dp p=a 



A 2 -a 

which, under the condition of no current flowing radially through the ring boundary, gives 
zero identically. Similarly, if we regard the anomalous Green function as playing the role 
of the Ginzburg-Landau order parameter in the theory of superconductivity, no-flow of 
supercurrent radially gives the average to be zero for . Furthermore, under the limit that 
the average of the product of two functions can be approximated by the product of two 
averaged functions, we obtain a set of reduced one-dimensional Gor'kov equations, with the 
coordinate variables being azimuthal angles 9 and 9' and the radii being set to the fixed 
value of the average radius. After appropriately renormalizing the Green functions by the 
inverse of the average radius (thus making the dimension consistent), we arrive at Eqs. Q. 
We remark that even if the ring has the geometry of a torus, the same averaged ID equations 
will result provided that the cross-section is relatively small, compared with the ring area 
and the coherence length. 



APPENDIX B: CRITICAL TEMPERATURE EQUATIONS: SUPPLEMENTARY 
DETAILS 

In this appendix, we supplement the details leading to the critical temperature equations 
in both clean limit and the disordered regime. 
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1. Clean limit 



Continuing from Eq. (Q, we make a shift in i: x — > x — m/2, under which the integral 
becomes 



E 

kez 



dx- 



p—i-wmk „i2nxk 



■(x - x ) 2 + n\[-iu ri 



(X + X ) 2 + fi\ 



(Bl) 



Jn 2MR 2 ^ i ^JL "^n 2MB? 

where x$ = — m/2. The integral over x can be performed using contour integration. For 
example, for k > and u n > 0, we can close the contour in the upper plane and evaluate 
the residues at x = x + \/2MR 2 (fi + iu n ) and x = —xq — ^2MR 2 (n-iuj n ). We end up 
with (noting the cancellation of the factor e~ l7Tmh ) 

7re i27rfc</>+i27rfc y/2MR 2 (ii+iuj n ) 



MR? 



(/I + iu n ) 



+ 



7re 



ix 



MR 2 



n 1 MR 2 

i2irkcf)-i2iTk^2MR 2 (fi-iiv n ) 



(jJL + iuj n ) 



MR 2 



(/! - lUJ n ) 



MR 2 + i3C V MR 2 

Including the contributions from the three other cases, we arrive at the self-consistency 
equation: 



(B2) 



\V\T 
2txR 



47T 



U>n>0 



MR 2 



(jJL + ioj % 



MR 2 



IXq 



MR 2 



(/I + iu ri 



+ E 



8vr e ^^2MR 2 {^ n ) cog 2^ 



(B3) 



where we have put the Debye frequency ujo as the upper cutoff in the Matsubara sum. This 
equation can be re-written as 



<2M 
/i 



V|TRe 



1 



w„>0 



1 + 



U n + 



i x 
AIR 2 



IXq 



(jj + lU n ) 



MR 2 



fc=i 



2e i2nk^/2MR 2 (n+iu Jn ) cog 3^ 



(B4) 



Furthermore, for typical temperatures we have + iuj n ~ yT^l + ^n/2/i). The second 
term causes the contribution of the fc ^ terms to be exponentially small for large R, 
i.e. there is a factor 

e -2nkR^2Mfj,w n /2fi 
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2. Large-radius limit 



For sufficiently large R, we can ignore the k ^ correction terms, in which case we have 



'2M 



\V\TReJ2 



1 



, n >o | Jl+i^ 



(B6) 



The third term in the denominator is usually much larger than the second term, as the 



chemical potential \x is much larger than the level spacing 1/2MR 2 , i.e., ^/2MR 2 /2 > 1. 
For typical values of lod and fi, we have up/ 1 H <C 1. Hence, we have 



'2M 



WD 



V\TRe 



1 



2u 



>2M\V\ 
\X 2tt 



U> n >0 tO n — IXQy J£r2 
lud /2-7rT 

Re Y. 



1 



1 \ _ ix-p / 2^ 



(B7) 
(B8) 



2ttT V A/fl 2 " 



The solution for T to this equation gives the critical temperature T c . Denoting by T c ° the 
critical temperature in the absence of flux (so that xo = 0), we have the corresponding 
equation 

n + 



'2M\V\ 



fi 2tt 



(B9) 



n=0 '" ' 2 

Taking the difference of equations flB7p and flB9p . we have 



w d /2itT c {4,) 

= Re 



1 



w D /27rT c ° 



„=o (n + D- ^yyil 



T — 

n=0 ^ 



(BIO) 



If we extend both upper limits to infinity, we should compensate by the difference (assuming 
uj d /T > 1), i.e., 

(w D /27rT°)+l 



(w D /27rT c (</>)) 



hi 



n 



(Bll) 



+i 



Therefore, we arrive at 



n.=o V (n + I) 



2ttT c O) Y Mil 2 



2^ n + 



In 



T c (0) 



(B12) 
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In terms of the digamma function ip(x) in Eq. (TTTj) . we have an implicit formula for T c (0): 



T c ° v \2J ^2 2vrT c (0) V MR 2 ' 

Using the quantities /i ~ Vp/2M, £ = vf/tiA , and p ~ -R/£o ; an d defining t((f>) = T c (0)/T c °, 
we can re-write Eq. (IB13[) as 

kf w =* G)- r ^ G- < bi4 » 

3. Finiteness correction 

What is the correction to T c (0) due to the finiteness of the radius? To address this 
question we need to take into account the k 7^ corrections to Eq. flB4[) . If we take 
Ud/[J> <C 1 and V f 1 + ~ ^//^(l + ^n/2/u), the self-consistency equation (IB4I) can be 
approximated as 



'2M\V\ 
fi 2tx 



[l 2tx 



u D /2irT 






Re Y. 

n=0 




1 




ix rT c ° 

2pT 


uj d /2ttT 






Re £ 

71=0 




1 


(n+|) 


ix rT c ° 

2pT 



2eJ 2^ V ^^(M+^n) cos27r A;0 j (B15) 

fc=i 



2M|VL ^ 1 f ^^^-^(n+i) 



^{l + ^2e^ 2MKV - fc ^ (n ^ ) cos2 7 r^}. (B16) 
fc=i 

Taking the difference between this equation and the version corresponding to <fi = 0, and us- 
ing the trick for converting the cutoff at the Debye frequency in a logarithm (see Eqs. flBllj) - 
(IB13p . we arrive at 



M0) = ^(~) - Re^.(i-|^) - 4^ {cos (2nkVnm) e^RF^ 1, |e"^] 

^ fc=i ^ 

-Re{ e ^V^ e -^HF[i-^,l^-^,e-^] cos 2ttW(1-^£) ) 1 ,(B17) 

where the dependence of t on r.h.s. is suppressed. We have extended the upper limit of the 
sum over n to infinity for the exponentially decaying terms (which introduces a negligible 
small error), and we have used the formula 

E ZTT = 2e_6/2 HF [- + a, 1, - + a, e" fc ] / (l + 2a), (B18) 

n=0 * 

where HF[a, b, c, z] is the hypergeometric function 2-^1(0, b; c; z). 
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4. Disordered regime 



In this subsection, we calculate G u (n) and K w and obtain the r.h.s. of the self-consistency 
equation (jHTj) . For simplicity, we assume that the potential u(r) is short-ranged, so that its 
Fourier transform u(q) can be treated as a constant, essentially independent of momentum 
transfer q. 

The one-particle self-energy can be obtained from summing one-particle-irreducible dia- 
grams (and ignoring the crossed diagrams), and thus we obtain 



S M = E w'"' 2 ^,. « -^gn(ul (B19) 

the real part of the self-energy has been ignored and we shall call t(u) the frequency- 
dependent scattering time. Under the condition that the Debye frequency is much smaller 
than the chemical potential, i.e. uj£,/ ' \x <C 1, and hence for the range of tu's that are relevant 
to superconductivity (i.e. < cod), we obtain 
1 1 



2r(u) 2t L 

where 



1 + 2 cos(2ttA;0) cos[27rfc ^2MR 2 /j] e ~ 2nk vW^ , (B20) 

fc>0 



27™ imp |u| 2 



(B21) 



2r y/2fi/MR 2 ' 

This form for t(u) is an approximant because it is calculated with G° rather than G. To 
improve the approximation we then use the exact disordered Green function G to calculate 



the same self-energy diagrams again (see, e.g. Ref. [23(), thus obtaining a self-consistency 
condition for t(u) (that can be solved iteratively): 



1 



2t(u) 2t q 



1 + 

fc>0 



2 cos(27rA;0) cos[2ttA; ^2MR 2 h] e - 27vk yfW^ ( 1+ ^Tr) 1 . (B22) 



As the k > terms are exponentially small, we can approximate r{uS) on the right-hand side 
by t to arrive at 



1 



1 

fc>0 



Y 2 cos(2vr£;0) cos [2tt k v / 2MR 2 fj] e'^y/W^ ( 1+ ^h ) 1 . (B23) 



2t{u) 2t q 

Returning to the kernel we re- write Eq. (I3T1) . following Gor'kov as 

K^m) = G u (m) G^{n 2 ) [l + C^m + n 2 )] , (B24) 
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where we have conveniently dropped the sign in G, and 

C u (ni + n 2 ) = n imp ^|u(g)| ! X(ni-g,n 2 + g) (B25) 

Q 

and the argument ri\ + n 2 in £ reflects the conservation of total (azimuthal) momentum of 
the two incoming (and outgoing) electrons in the ladder diagrams after disorder averaging. 
For the purpose of evaluating self-consistency equation ( |32|) . our goal is to obtain for 
n 2 = —rii + m. By eliminating K from Eqs. flB24j) and (IB25I) by we have an equation for 
r ■ 



which gives 



where 



C u (m) =n imp J2H 2 G UJ (n' 1 )G^{-n , 1 + m)(l + C UJ {m)), (B26) 



£ »« = ott< (B27) 



A u {m) = n imp \u\ 2 G u {n[) G^(-n[ + m). (B28) 

The self-consistency equation (I3"2"j) then becomes 

1 _ V 1 AJjn) _ V 47tto ^ AJjn) 

~ (3L^n imp \u\ 2 l-A w (m) /3L y/2fi/MR 2 ^ 1 - A w (m) ' 

where we have used Eq. (1B21I) in the second equality. As with the evaluation of Eq. (1B1I) in 
the clean limit (see Appendix [B]) , we obtain as 



4° _i_ A k>0 



—Re 
2r 



, V cos(27rkJ2MR 2 ij,)Re — , (B30) 

11/ u fc>Q 11/ u 



where A = ((f) — m/2)y jj^z, r\ = (1 + 2t|u;|)/(2t|u;|), and r is a shorthand for t(o>). 
To illustrate the corrections arising from the finiteness of the radius, we re-write A^ as 

A _ (l + 2r|^|)(l + qJ-2r^A 
w " (1 + 2t|o;|)2+(2tXo) 2 ' {B6l} 

where a w (^C 1) and b w (<£. 1) are 



^ 2e ~ 2nk ^W^ r i cos(27rk^2MR 2 fi) cos(27r£;0), (B32a) 

fc>0 



2e^ k y^^ T ' cos(27rk^2MR 2 fi) sin(2vrfc0). (B32b) 

fc>0 
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As the correction to l/r(u) is exponentially small [see Eq. (1B23I) ]. in the exponentials rj can 
be approximated by 

= 1 + s^f- < B33 > 

On the other hand, for 1/t(u) not in the exponentials we shall approximate it by retaining 
the leading correction [see Eqs. ( 1B23I) and ( 1B32al) ] via 

^-(1 + a u ). (B34) 



2r(u) 2t 
The quotient A u /(1 — AJ) then becomes 

a 1 _|_ \a — 2r^X 1 I , r _ 2r^X 1 



(B35) 



1 " 4- 2rM + TT ^(2rX )2 - [a. - fg-g] 2r| W | + j^firX^ 

We remark that from this equation one can obtain the equation for the critical temperature 
for arbitrary mean- free path l e . 



a. Large-radius limit and arbitrary disorder 



To illustrate the remark made above, we consider large-radius limit in which we ignore 
corrections due to the finite radius. In this limit Eq. (1B35|) becomes 

A„ 1 



1 — A^ 2tq\uj\ 
Substituting this into Eq. (1B29I) we have 
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In the clean limit tqT c ^> 1, we see that this reduces to Eq. ( 1B8I) . which we obtained in the 
absence of disorder. But here Tq is arbitrary. Subtracting Eq. (IB37j) from the corresponding 
equation at = (i.e. X = 0) and using the trick to get logarithm of the ratio of the 
critical temperatures as we did in Eqs. (1Bllj) - (1B13l) . we have 
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By making the partial fractions of the first term in the summation and by using the definition 
of the digamma function fTTTj) . we arrive at 



In 



T c ((f>) 
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in 



a 2 — x 2 L 



— a + \J a 2 — x 2 ( 1 + ot + a/« 2 — x 2 



a + a/" 2 - X 2 , / 1 
^ V> 



a 2 - x 2 



(B39) 



where, for convenience, we have defined a = l/47rroT c (0) and x = Ao/7rT c (0). This is the 
equation for the critical temperature at arbitrary disorder for large radii of rings. 



b. Finite-radius corrections and strong disordered regime 

Now we return to the corrections due to the finiteness of the radius. In the strong 
disordered limit (i.e. r u <C 1), Eq. (1B35I) becomes 

1 - 2rb^X Q 



l-Au 2t[\u\ + 2tX 2 ] 2t[\u\ +2tX 2 ] 
This is further approximated, using Eq. ( 1B34I) as 



(B40) 



1 2a w -2r 6 w X , D .,> 

+ n — r i i i o — ^T' y BA1 > 



1 - A w 2r [\u\ + 2t X 2 ] 2t [\u\ + 2r X 2 ] 

where we have used Eq. (1B34|) . The first term on the r.h.s. represents the bulk term and is 
the Little-Parks term in the dirty limit. The second term takes into account the finite radius, 
and contains flux dependence in period h/e. The equation for the critical temperature will 
then contain the digamma function and hypergeometric function with real arguments, in 
contrast with the clean limit. This makes the transition to normal state continuous, as 
the solution for the critical temperature is unique and the assumption of vanishing order 
parameter used in the linearized self-consistency condition is valid. Therefore, there can be 
quantum phase transitions tuned by flux and/or radius. 

As ro|Xo| <C 1 in the strong-disorder limit, the ratio A^/il — A u ) becomes (by ignoring 
the b term) 

^ 1 | (B42) 

1 - A„ ~ 2r [M + 2t X 2 ] 2t [\u\ + 2t X 2 ] ' 

From this we can obtain an equation for the critical temperature T c in the presence of flux 

$, as we have done in the clean limit (see Appendix IB 3h By inserting Eq. (1B42|) into the 
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self-consistency equation (1B29j) we obtain 



D 
2nT c 
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/i 27T 



n=0 



Subtracting from this the corresponding = equation, i.e., 
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we obtain the following implicit equation for the critical temperature: 
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where we have used the fact that when n « lje>/2iiT® or higher, + | + x) ~ 1/n, the 
sum of the series being approximated is a logarithm, and we have extended the upper limit 
of n for the correction terms to infinity. By using formula (1B18I) for the hypergeometric 
function, we finally arrive at 
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where we recall that X = (0 — m/2)^J~^^. If we use /i = Mvp/2, l e = vpTo, and 
£0 = vf/ttAq (with A = rT c °), we have 

lnt = Ki) ~ ^2 + " V) 2 ) + E 4e ^ cos(2vrA; v ^M^) 
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We can re-write this equation using the relations /i = Mv F /2, A = Vf/k£o and l e = VfTo, 
and thus arrive at 
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Finite-radius and arbitrary disorder 



It is also possible to include corrections from the finiteness of the radius, and derive 
the equation for the critical temperature, as we did in the previous section. The results 
will contain the r.h.s. of Eq. (IB39I) as the zeroth order, as well as corrections due to the 
finiteness of the radius in terms of hypergeometric functions. However, the formulas are too 
cumbersome and we do not list them here. 



APPENDIX C: HEURISTIC ARGUMENT TO THE EMERGENCE OF h/e PE- 
RIOD OSCILLATIONS 

1. Cooper problem on a ring 

Let us consider Cooper's problem on a ring with a flux $ threading through it. The 
orbital wavefunction can be written in the same form as in Eq. 03]). Assuming that ip(x, x) ~ 
e i2nmx/L^ an( j ^ e e i ec tron-electron interaction is factorizable, i.e. V n>n i = AU n U*,, we 
can write down the time-independent Schrodinger equation in terms of the wavefunction 
amplitude 

2e n (<p)a n + ^ V n ^a nl = Ea n , (CI) 

\n'\>np 

where 

26n(0) = K n + <P) 2 + (m + n- <pf] , (C2) 

and hence we arrive at 



A ^E-2e n {<PY 
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The integer m has to be chosen to minimize the ground state energy. Thus, we see that the 
solution is periodic in <p with period 1 (or in $ with period h/e, namely the single particle 
flux quantum) and that the Little- Parks period h/2e is exact only in the large- R limit, in 
which case the summation over n can be replaced by an integral and hence the period of 4> 
is 1/2. This illustrates the role of the single particle oscillation. However, the amplitudes 
of the single-particle and Little-Parks oscillations have to be evaluated within a microscopic 
theory, as has been done in the main text. 



2. Flux oscillation: an instanton approach 

The instanton picture provides us another, heuristic, view of the emergence of /i/e-period 
oscillations. The argument presented here is intended just to give some intuition. We refer 



to Rajaraman [27| for a pedagogical review of instanton techniques. A single instanton 

tunneling from through a potential to 2ir (which is identified as 0) on a circle gives an 
amplitude 

lim (27r|e-^ r |0) (1 o) = e" 50 Jir™ 1/2 e" wr/2 , (C4) 

where So is the classical Euclidean action (in the absence of any flux threading through the 
circle), J is a Jacobian factor, K is a constant independent of r as r — > oo, and u> is the 
harmonic oscillator frequency near the bottom of the trapping potential. Now, if there is a 
flux $ threading the ring, the amplitude will also acquire an additional phase factor e l2lT ^, 
where = $/$ an d $o = h/e, if the particle carries charge e. For a charge 2e particle, the 
phase factor would be e 1471 "^. 

We consider a ring size that is large enough that the electrons in a Cooper pair generally 
traverse the ring together, and rarely split so as to wind separately around the ring. The 
ground state energy will give a qualitative estimate of the critical temperature of the asso- 
ciated superconductor. The total contribution of the amplitude for the instanton associated 
with the Cooper pair is then 

lim (27r|e- ffT |0) = u^e'^ 2 Y -J—UKre'^V^e^-^, (C5) 

711,712 

which gives 

lim (27r|e- Hr |0) = cu 1/2 e^ T/2 exp (2 JKre^ cos4tt0 N ) ~ lim (27r|£ )(£o|0}e" EQr . ( C6 ) 
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This gives a ground state energy of Eq = ^ — 2 JKe" cos 47r</>, which reveals the period 
/i/2e, the Little-Parks period in critical temperature. The fact that this gives a dependence 
on flux being sinusoidal rather than quadratic results from the lack of accounting of the 
other many-body electrons and the existence of a condensate. 

Now, occassionally (in the sense of contributing Feynman paths) the electrons in a Cooper 
pair separate and circumnavigate the ring (relative to one another) before re-associating. 
The contribution of such processes to the ground-state energy can be estimated via the 
instantons and anti-instantons of such events associated with them. A single instanton 
involving one electron going from to 2tt has the amplitude 

lim (27r|e- HT |0)no) = -e~ s ^ J e K e ru l J 2 e-^ T/2 ^ , (C7) 

where the subscript e indicates a single electron rather than a Cooper pair, and the additional 
minus sign comes from the exchange of the two electrons. Summing all the instanton and 
anti-instanton processes, we arrive at the contribution to the ground-state energy from the 
two electrons 

E . e = u e + AJ e K e e- So - e cos 2vr0, (C8) 

which results in the emergence of an h/e contribution to the period of oscillation to critical 
temperature. We remark that it is owing to the separation of the lengthscales and hence 
time scales that we can separate the Cooper-pair and single-electron contributions. The 
amplitudes of the two oscillations are related to the respective actions, and for large radius, 
the amplitude corresponding to single-particle oscillation is expected to be small, due to the 
binding resulting from the attractive interparticle interaction. 
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